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Superstatistics describes statistical systems that behave like super- 
positions of different inverse temperatures /3, so that the probability 
distribution is p(ei) oc /(/3)e^'''''d/3, where the 'kernel' is 
nonnegative and normalized (//(/9)d/3 = 1). We discuss the re- 
lation between this distribution and the generalized entropic form 
5 = s{pi)- The first three Shannon-Khinchin axioms are as- 
sumed to hold. It then turns out that for a given distribution there 
are two different ways to construct the entropy. One approach uses 
escort probabilities and the other does not; the question of which 
to use must be decided empirically. The two approaches are re- 
lated by a duality. The thermodynamic properties of the system can 
be quite different for the two approaches. In that connection we 
present the transformation laws for the superstatistical distributions 
under macroscopic state changes. The transformation group is the 
Euclidean group in one dimension. 

thermodynamics | entropy | classical statistical mechanics | correlated systems 

Superstatistics has been introduced as a way of sys- 
tematically handling statistical systems that can be seen 
as superpositions of various Boltzmann distributions. One 
may think of a superstatistical system on a microscopic and 
a macroscopic scale. At the microscopic scale the system re- 
laxes towards thermodynamic equilibrium and follows a single 
Boltzmann statistic with a well defined local inverse tempera- 
ture /3, i.e. the local probability of finding the system at some 
energy e is proportional to exp(— /3e). On the macroscopic 
scale the local inverse temperatures fluctuate, with the fluc- 
tuations governed by /(/3) Q- Following [1], the generalized 
superstatistical Boltzmann factor reads 



B{e) 



dl3f{p)e 



-fie 



where is called the superstatistical kernel. The kernel is 
non-negative and normalized, i.e. dp /(/?) = 1. The prob- 
ability of finding the system in one of W discrete (for example 
energy) states, e = {ei}J!^i, is 



Pi = -^Biei) 



[2] 



Z being the superstatistical generalization of the partition 
function. The expectation value 



U 



E 



Piti 



[3] 



is the internal energy of the system. 

Superstatistics has found many applications ranging from 
hydrodynamic turbulence [H O |4j , complex networks [5] , and 
pattern formation [6] to finance [71 18] . Besides these practical 
aspects, superstatistics has the potential for providing a struc- 
tural foundation for non-Boltzmann statistical mechanics. 

Different methods have been proposed for establishing the 
mathematical methodology of a kind of statistical mechanics 
compatible with superstatistics [51 IIUI [TTl 1121 [T5] . The meth- 
ods that use generalized entropic forms are the Tsallis Souza 



(TS) approach [9] and one described in [121 [13] which we refer 
to as the HT approach. Both use a maximum entropy prin- 
ciple (MEP) to reconstruct generalized entropies having the 
structure 

w 

S[p]=^s(pO [4] 

!=1 

from given superstatistical distribution functions. Here p 
stands for the set of probabilities {pijjli. The sum struc- 
ture is chosen in analogy to Boltzmann-Gibbs (BG) entropy, 
where sbg{x) — — a;log(a::). The function s is continuous and 
concave and has s(0) = 0. These properties correspond to the 
first three Shannon-Khinchin axiom^ [141 115j . In this paper 
we address the question of whether the TS approach is the 
unique way to reconstruct entropies for superstatistical sys- 
tems. The answer is no - there exist possible alternatives. In 
fact, if we did not have the first Khinchin axiom we would 
have infinitely many possibilities. We show however, that as- 
suming this axiom we have only two possibilities remaining. 
The axiom requires that the entropy be a continuous func- 
tion of the probabilities p only (compare footnote [5]| . It turns 
out that one possibility corresponds to the TS [9] and the 
other to the HT approach [HI [13]. We show that the HT 
case cannot be seen as a simple limit of the TS case. For 
a large class of distribution functions such a limit yields the 
BG case, not HT in general. The thermodynamic properties 
of the entropies corresponding to the two approaches will in 
general be quite different. Which approach is the correct one 
to use in each case has to be decided empirically. Superstatis- 
tics - as originally presented- does not make any assumptions 
about how the superstatistical Boltzmann factor B or the ker- 
nel / transform under thermodynamic state changes. On the 
macroscopic scale, a superstatistical system can be described 
by a macro-state a, in complete analogy to a classical ther- 
modynamic system. These states are characterized by a set of 
macro- variables, e.g. internal energy, volume, temperatures, 
a = {U,V,T, ...}. We present the transformation group of 
the superstatistical kernel / under changes from one state to 
another, a — !■ a'. 



Entropy reconstruction from distribution functions 

The usual procedure when using the MEP is to maximize 
a given entropy under given constraints, represented by La- 
grange multipliers: 

^[p]^S[p]~a(^p,^l^ ~p(^p,e,-U^ . [5] 

Given s{pi) and [/, this yields the distribution function pi, 
along with the values of the Lagrange multipliers a and 13. 



"'^ Superstatistical systems are usually not in equilibrium but in stationary non-equilibrium. E.g., 
tW\nk of a gas or liquid between two plates at different temperatures Ti and T2 . 
^Sliannon-Khinchin axioms: (i) Entropy is a continuous function of the probabilities pj only, i.e. 
s should not explicitly depend on any other parameters, (ii) Entropy is maximal for the equi- 
dlstribution pi — 1/W. - From this the concavity of s follows, (iii) Adding to a system a state 
numbered W + 1 with p\Yj^i — does not change the entropy of the system. - From this 
s(0) ^ follows. 
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In the case of entropy reconstruction the distribution func- 
tion of a superstatistical system is given. The system is in 
a specific physical stat43 o"- The various quantities such as 
Pi , the internal energy etc. are given and are labeled 

by ((j). Since p'"^' is a superstatistical distribution, we have 

pf'' = _B(e*°'')/Z''^'. The quantity s{j>i) and the entropy can 
now be computed for the state a. 

Entropy reconstruction - the HT approach. We begin with 
the simpler case and demonstrate the entropy reconstruction 
using the HT approach [121 113) . which starts by maximizing 
the functional of Eq. ((5]). At this point we use the available 
information about the specific state a. The condition 







[61 



leads to 



s'ht{p. 



(a) . n((T) (cr) 



[7] 

Since i3''^'(e) is a monotonic function of e the inverse function 
L defined by 

exists. Note that for the specific state cr, Z^"^^ is just a pos- 
itive real number. Otherwise the inverse would be hard to 
calculate. It follows immediately that 



[91 



for all j. Equations ^ and Q yield the differential equation 

s'ht(p1^') . [10] 

This allows us to reconstruct the entropic form by integration: 

aHT{x) = a^-'^x + p^''^ r dyL{y) . [11] 
Jo 

We have now derived an entropic forirQ sht which, when 
used in the MEP, Eq. (O, reproduces p^^\ a^'^\ and P^"'^ 
as a solution, given the state characterized by C/'-'^-' and ef\ 
As a result of the first Khinchin axiom this entropic form 
sht characterizes the superstatistical system in general, i.e. 
not only for the state a but also for the other possible states 
cr', corresponding to different values of f/''' ' or ef This 
has consequences for the thermodynamical treatment of the 
system. In particular, transformation laws for the supersta- 
tistical Boltzmann factor B and the kernel / under changes 
of the state cr can be derived, as will be shown below. 

Entropy reconstruction with escort distributions - the TS 
approach. We now generalize the above idea by modifying the 
energy constraint. We compute a different 'energy' U* by us- 
ing a different set of 'probabilities'. These probabilities can 
be written in the general form 



u{pi 



J2j -^(pj) 



[12] 



The quantity Pi is often called escort probability. The escort 
energy U* is identical with the internal energy U if and only 
if u{pi) = Pi. We do not attempt to provide a physical in- 
terpretation of the escort probabilities. We start with the 
functional 



<f [p] = S[p] 



E 



Pi 



■ U' 



[131 



For taking the derivative of P, we define the following notation 
using a function Q of two variables: 



u{p,) 

Ul{pi) 
U2{pi) 



^l-X^y] \x—pi-^y — s(pi) 



a_ 



[14] 



=s(Pi) 



Maximizing the functional in Eq. (|13p we have d^/dpi — 0, 
leading to 

s'{p^)=a + |3''{e^~U'^)[ul{p^)+U2{p^)s'{p^)] , [15] 

where /3* is the rescaled parameter P/ '^^u{pi). Solving for 
s' , we get 

cv + /3*(e, - (7*)ui(pO 



s'(pi) 



- P'{e,-U*)u2{pi) 



16 



As in the HT approach, we consider from this point on a 
specific state a of the superstatistical system and use the spe- 
cific values 0!*°^^ for reconstructing the associated 
entropy s. However, for the sake of readability we suppress 
(cr) whenever that cannot cause confusion. Using Eq. ((9} as 
above, we insert — L{pi) into Eq. (|16p and get 

a + P'iLix) - U*)ui{x) 



s'{x) 



[17] 



I- P*{L{x)-U'-)u2{x) ■ 

This differential equation in s can now be solved, so that 
for any suitable function u an entropic form s is obtained 
which reproduces the superstatistical distribution function 
p(<T) ^ under the MEP. 

It is easy to see that such entropies could explicitly de- 
pend on U* , except for the first Khinchin axiom, which rules 
out such dependence. This reduces dramatically the possibil- 
ities for choosing it. In fact, the first axiom can hold only 
if both ui{x) = ri and U2{x) — r2 are constants, implying 
^(Pi) = '"iPi + T^2s{pi)- This now allows us to absorb U* into 
the Lagrange multipliers. By defining new constants 

, a-P*U'ri J ,3 P'r-, 



l + P*U*r2 
we can simplify Eq. (|17l) to 

s'{x) = 



and P = 



1 + P*U*r2 



a + PL{x) 
1 - PuLix) 
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[191 



with v = r2/ri. This equation can now be solved by simple in- 
tegration. Note that since Pi remains invariant if both ri and 
r2 are rescaled by the same factor we can use ri = 1 without 
loss of generality. With u{pi) = Pi + vs{pi), the TS approach 
is recovered. Note that more general constraints than those in 
Eqs. ([5]) and (|13p are ruled out by the first Khinchin axiom. 

The two entropies 

We see that only two entropies remain, one obtained with 
ordinary constraints (HT approach), and the other by em- 
ploying escort distributions in the functional of Eq. (|13|) (TS 
approach) . 

• HT entropy: The approach without escort distributions 
results in Eq. (fTTj) 



[20] 



sht{x) = P I dyL{y) + ax 
Jo 

which yields the HT entropy, [HlfTH] . 



■ The state cr does not change during the reconstruction. 

In rT2]fT3l A(x) — —a — fiL{x) is called the generalized logarithm. 



TS entropy: The only possible entropy involving escort 
distributions results from Eq. (|19|) . which gives 



sts{x) = / dy- 
Jo 1 



the TS result for the entropy, and Pi 



(Pi + 



vsrsiPi)) / ^j{Pj + i^^TsiPj)) for the constraint. In [3] 
j3 is set equal to one. 

Note that the TS approach has the freedom in choosing v 
(each choice of v leading to a different entropy Sts), whereas 
in the HT approach the entropy is fully determined by the 
distribution function alone. The question arises of whether 
there exists a way to unambiguously determine the value of 
V. That is addressed in the following section. 



Duality of the two approaches 

Using Eqs. H20p and (|2ip we have for the functional relation 
between the two entropies 



i (q/3 - a/3) + s'ht{x) 



[221 



Suppose we reconstruct an entropic form a from state a. We 
now use this form in the MEP for another state a' . This could 
for example be a state with a different internal energy. For 
both approaches, the HT and the TS, solving the MEP for 
state a' leads to distribution functions of the form 



where for the HT approach we have 



and for the TS approach 



and h - 



[231 



[241 



and 6=£^i±i^ . [25] 



From Eqs. (|24|) and (|25p and the requirements lim^-^o a = a 
and lim„_>o /3 = /3 we get 







ajS — a/3 and 1 = ^ ~\- va 



so that Eq. (|22p simplifies to the relation 
1 1 



[261 



[27] 



Remarkably a and /? drop out. This relation can be expressed 
in terms of so-called generalized logarithms (g-logarithms) for 
the HT and the TS approaches. Generalized logarithms A{x) 
have been widely used in the context of generalized entropies, 
e.g. |16l I17| . A is an increasing monotonic function with 
A(l) = and A'(l) = 1. Equation ({27} suggests that given 
s'j-g{x) = for some x = 2:0, ShtI^^o) is also 0, and further 
Sts{xo) — s'lij'{xo). If we define a constant c = —xo s'^rj^l^xo), 
the g-logarithms associated with the generalized entropies can 
be written as 

1Vht{x) ^ --s'„t{xox) , kTs{x) = --s'ts{xox) . [28] 
c c 



With these g-logarithms Eq. (|27p becomes the duality relation 
Aht{x) = Ats{x) , [29] 



[21] where the duality operation * is defined by 



A* = 



+ cu 



[301 



Note that Eq. (|30l) possesses a symmetry: the equation 
is invariant under interchanging A (A = Aht) with A* 
(A* = Ats) and simultaneously changing the sign of v. Ob- 
viously, applying * a second time (exchanging Art •<-> Ats 
and f <-> —v) yields the identity, as has to be the case for any 
dualitjQ. In Eq. O we see that the unique value of v, where 
A* and A and hence sht and sts can coincide, is v = 0. On 
one hand, in the v ^ limit the HT and the TS entropies 
are obviously identical. On the other hand superstatistical 
systems are generalizations of BG statistics. It is therefore 
natural for the two approaches to coincide for the BG case. 

Example: duality and power laws. For a large class of dis- 
tribution functions the duality relation of the associated g- 
logarithms includes the condition 



A*(x) = -A{l/x) 



[311 



given in [TS]. More precisely, a large class of families of g- 
logarithms is closed and the usual logarithm is self-dua0 un- 
der this map. By requiring that both conditions, Eq. (|30p 
and Eq. ((311), hold, i.e. by inserting Eq. ((ST)) into Eq. ((SO)), it 
is possible to derive explicitly a generic form of g-logarithmfl 



A{x) = 



[321 



log(a;)) 



with A G (0, 1] and h{x) a monotonically increasing func- 
tion h : [-cx),cxj] [-1,1], with h{0) = 0, h'{0) = 1, 
h{x) = —h{—x), and limx^oo h{x) = 1. Each pair h and 
A defines a family of g-logarithms parametrized by v; each 
such family A,^ has the properties (i) Al = A_^ and (ii) 
lim^-t-o Ay = log. The first property states that the g- 
logarithms associated with the TS and HT entropies are dual 
to each other. The second property states that in the i/ — > 
limit all these families reproduce the BG cas43- This result 
is quite general: if sht is associated with a g-logarithm A^i 
of the form Eq. (|32p for some value u' then sts is associ- 
ated with A_y/, with v' fixing the value of f in the escort 
probability. 

The specific choice of h{x) — tanh(a-) and A = 1 yields the 
family of g-logarithms 



A{x) 



■ tanh(^ log(x)^ 2 



log„(2;) 



[331 



where \og^{x) = (a:^"^ — 1)/(1 — q) and u = 1 — q. We see 
in a concrete example how u is uniquely determined by the 
parametrization (here q) of the family of given distribution 
functions. It is remarkable how requiring both, Eq. H3Up and 
Eq. (|31|l . automatically reproduces the condition necessary to 



^ The duality * acts on pairs {A, v) in such a way that (A, iv) * — (A* , i^* ) with A* as in Eq. 
f30] and /y* ^ -ly. All elements (A, 0) are seif-dual, i.e. (A, 0)* = {A, 0). This allows us 
to explore the HT-TS duality in many different directions that go beyond the scope of this paper. 
"^Self-duality here means log*{cc) — log{a:). 

''^Without loss of generality we can set c — 1 since c trivially rescales the escort parameter v. 
^ Eq. |31| as the starting point for producing *-closed families of g-logarithms with properties (i) 
and (ii) is not the most general Ansatz. Other possibilities will be discussed elsewhere. 
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recover Tsallis entropy [18] within the TS-approact0. More- 
over, the duahty = written in terms of g-logarithms 
is just the well-known duality log* = log2_q. The associated 
HT and TS entropies are 



q-l 

ShtIp] = -EJ<f'rfx-log,(-^ 



1-9 
9(2-9) 



with 2:0 = Q^''*^ defined as above by s'j'g{x)\^^^^ = 0. 
Clearly Sts gives the Tsallis entropy [19] while Sht corre- 
sponds to the entropy for power-laws discussed in |12l 113) . 



State changes - the superstatistical transformation 
group 

As we have seen before, Eq. (|23p determines how the su- 
perstatistical Boltzmann factor B and the kernel / transform 
under state changes ct — a'. These transformations are 



^Z^"' and B<"')(e) = ^B^"'(a-f 6e) 



[351 



where 2; is a ratio of two normalization constants for a and a' 
respectively. For the transformation law of the superstatisti- 
cal distribution function we find 



-^(a + b€) 



i/J-d/3/(-)(/?) 



[36] 

where we substituted /?' = bp in the second line and used Eq. 
([35} in the third line, together with Eq. ([l}. By comparison, 
/ transforms according to 



Z 



1 1 



P' 



[371 



The value of z is fixed by the normalization condition for the 
kernel, J^dPf'^"'\p) = l. 

For describing the transformation group we introduce the 
following notation: state a is characterized by a = and 
6=1. For state cr', a and 6 take other values. We now 
indicate the parameters (a, 6) in the kernel through 



j(<t) _ ^(0,1) 



Equation ((37|) can now be written as 



[381 



[391 



From Eq. (|36|) we see that under the state change a ^ a' the 
energy e undergoes the affine transformation e ^ a + be. We 
define three operators $1, $2, and 11 representing translation 
of e, dilatation of e, and normalization of /, respectively 



$i(a)/(/3) = e-''"/(/3) 
*2(6)./(/3) = 

n/(/3) = j^WaW) 



[401 



$1 and "1?2 form groups 

$2(6)$2(b') 



[41] 



$1 (a + a') 
^2 {bb') , 

with the identity elements $1(0), and $2(1), and the inverse 
elements 



<E>i:'(a) = *i (- 



■1(6) = $2 (I 



[34] and $2(6) do not commute: 



$2(6)-l'i(a) = -l'i (^j*2(6) 



[421 



[431 



while the projection operator 11 commutes with $2 • We finally 
define the group 

G(a, 6) = n $2(6) $i(a) , [44] 
with the following group-composition rule and inverse element 
G{a,b)G{a',b') = G(a' 4- a6', 66') 
G-i(a,6)=G(-^,i) . 

Equation (|39p can now be written as 

/("•'"(/3) = G(a,6) /("•i'(/3) , [46] 
or more generally, 



[47] 



The generators of $1 and "1>2 are given by gi — — /3 and 
(?2 = —(1-1-/3^). We have a representation of the Euclidean 
group in one dimension. 

These transformations apply to infinitely many different 
/'s. For example, if /'■"'i' leads to a g-exponential distribution 
with a specific value of q, under these transformations all pos- 
sible will give g-exponentials with the same value of q. 
The functional form of distribution functions is preserved un- 
der the transformation group. More generally, for every type 
of superstatistical distribution function, we can identify a rep- 
resentative /'"'i' and call it the 'progenitor' function. From 
this a sequence H(/) of can be generated by applying 
the transformation group. For every type of function / there 
exists a unique sequence to which it belongs. These sequences 
define equivalence classes for the superstatistical kernels. Sup- 
pose experiments on a superstatistical system produce a set 
of functions Q — {/o-^jfLi for different states ai, arbitrar- 
ily chosen. If f2 C H(/) for some fixed sequence S(/), then 
state changes of the system transform consistently under the 
transformation group Eq. (|44|) and an entropy S = ^(Pi) 
exists such that s is invariant under state changes of the su- 
perstatistical system. 

Discussion 

By studying the possibilities of reconstructing entropies from 
superstatistical distribution functions we found that there ex- 
ist only two distinct frameworks, the HT approach (using only 
normal constraints) and the TS approach (employing escort 
constraints). A duality between the two approaches is shown 
to exist. This indicates that superstatistical systems can be 
explored in two ways which are dual to each other. Under 
certain circumstances the duality allows us to determine the 
possibilities of functional forms of generalized logarithms and 
thus of entropic forms, as explicitly shown for systems with 



If a q-exponential is observed in the superstatistical system for some particular value of q, then 
the entropy reconstruction in the TS-approach produces different entropies 57^5 depending on the 
choice of h' in the TS-MEP. Each of these entropies can be used to recover the povyer-law (q- 
exponential distribution) one has used to construct the entropy. Yet only for the particular choice 
u — 1 — q does one get ^TS iPi) Tsallis entropy Sq [p] for the particular case [9]. 



characteristic power laws. For systems governed by a par- 
ticular distribution function it has to be decided empirically 
which of the two frameworks applies; these systems may dif- 
fer considerably in terms of their thermodynamic properties. 
We studied how the superstatistical kernel transforms under 
macroscopic state changes within the two approaches. We 
presented the respective transformation laws and identified 
the corresponding transformation group to be the Euclidean 



group in one dimension. This indicates the existence of a 
remarkably simple mathematical structure behind thermody- 
namical state changes of superstatistical systems governed by 
an invariant entropic form. 
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